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SUMMARY 


Differential equations of bending of an idealized H-section beam 
column were derived for a nonllnearly viscoelastic material whose 
mechanical properties are analogous to a model consisting of a linear 
spring in series with a nonlinear dashpot whose strain rate is propor- 
tional to a power of the applied stress. The resulting constant stress 
or load creep curve consists of a straight line, the slope of which can 
be considered as the secondary creep rate of a real material. 

The equations derived were used to obtain the creep-bending deflec- 
tions of a beam in pure bending and of a column with initial sinusoidal 
deviation from straightness. The results of the analysis of the simple 
beam showed that the deflections vary linearly with time. The analysis 
of the deflections of the column, accomplished with the assumption that 
the original shape of the structure was maintained at all times, showed 
the existence of a finite critical time at which the deflections become 
indefinitely large. The critical time decreases rapidly with increasing 
axial compression and column inaccuracy. 


INTRODUCTION 


The need for methods predicting the behavior of structural components 
at high temperatures is becoming increasingly urgent, particularly in the 
fields of aircraft structures and propulsion. During the past 40 years 
considerable attention has been paid to the fundamental constant-stress 
and constant- load tensile creep behavior of materials. However, with few 
exceptions, it is only recently that results of the investigation of the 
creep behavior of beams and columns have been presented (see, e.g., 
refs. 1 to 9)" It is the purpose of the present report to apply to the 
problem of the creep behavior of beams and columns a stress -strain-time 
relation which can be considered as a generalization of the relation 
obtained between the strain rate and stress for a Maxwell linearly visco- 
elastic model consisting of a spring connected in series with a dashpot 
(see fig. 1 and refs. 6 and 8 to 10 ). While the original Maxwell model 
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consists of a spring and dashpot, in each of which the strain is linearly 
related to the stress, the model considered herein consists of a linear 
spring coupled with a nonlinear dashpot whose strain rate is proportional 
to a power of the applied stress. The corresj>onding stress -strain-time 
relation can be adapted to the creep characteristics of structural 
materials, such as aluminum or steel, by assuming that the usual creep 
curve indicated by the dashed ciirve in figure 2 (a) can be replaced by 
the solid straight line also shown in the figure (see refs. 1 to 3 )- 
Hence, the idealized creep curve takes into account the actual secondary 
stage of creep, approximates the initial elastic or elastoplastic stage 
and the primary creep stage, and ignores the final stage. It is intended 
that the effective spring modulus, which defines the strain Intercept of 
the Idealized creep curve, together with parameters defining the action 
of the nonlinear dashpot, be determined from experiments at constant 
temperatures . 

In order to sln5)lify the calculations involved in the analysis of 
the creep behavior of beams and columns, the present investigation has 
been confined to the study of the behavior at constant temperatiare of an 
idealized H-beam, the cross section of which consists of two equal concen 
trated areas connected by a thin web of negligible bending resistance 
(fig. 3). Two differential equations for the determination of the deflec 
tions of a beam under combined axial and lateral loads were derived, one 
equation being applicable when the stresses in both flanges are in com- 
pression, the other when one flange is in compression and the other in 
tension. These equations were used to determine the deflection-time rela 
tions for a beam in pure bending and for a simply supported column whose 
axis before loading had the form of a half-sine wave. 

The author is Indebted to Professor H. J. Hoff for his guidance and 
criticism, to Mr. S. A. Patel for his assistance, and to the National 
Advisory Committee for Aeronautics for sponsorship of the research 
reported in this paper. 


SYMBOIS 


A total flange area of idealized H-section 

effective elastic modulus 

E- jE,. effective moduli for compression and tension, respectively 

XC Xu 

F = w/h 

Fc = w-cA 
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Fi = Vj^/h 






h 

I 

L 

M 

Mo 

m 


P 
t 
t. 


cr 


w 


Wi 



A 

total nondimens ional deviation of axis of loaded column from 
x-axis at t = 0 

amplitudes of Fj, and Fj^, respectively 

+ frp 

■^o 

i/[l - (ff/ag)] 

distance "between flanges of idealized H-section 

moment of inertia of idealized H-section, Ah^/k- 

length of beam 

bending moment 

constant bending moment 

exponent in viscosity term 

expionents for compression and tension, respectively 

axially compressive load 

time 

critical time 
deflection due to loads 

time -dependent deflection, accrued for t > 0 

initial deviation from straightness 

time -independent deflection due to loads at t = 0 

total deviation of axis of loaded column from x-axls at 
t = 0 


X 


axial coordinate of beam 
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z = 

6 strain 

strains on concave and convex sides of beam, respectively 
(positive in compression) 

strain at t = 0, 

X viscosity coefficient 

viscosity coefficients for compression and tension, 
respectively 

p radius of curvature 

cr stress 

stresses on concave and convex sides of beam, respectively 

°-K = rt^E^l/AL^ 

cTq constant stress 

a average compressive stress, P/A 

T = |Ei(2a)7[U(cTE - ■^)]} t 

time parameter corresponding to z = 1/2 
■^cr critical value of time parameter . 

n=a( )/st 

C )x = 

EQERIVATIOIT OF DIFFERENTIAL EQUATIONS OF BENDING OF 
AN IDEALIZED H-SECTION 
Stress -Strain-Time Law 

The relation between the strain rate and the stress corresponding to 
the model previously described is 
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5 = * (o“A) 


( 1 ) 


in which e and a, respectively, are the strain and stress, is 

the effective modulus, m and X are parameters defining the viscous 
behavior of the material, and the dot over a symbol indicates differenti- 
ation with respect to time t. The three parameters E]^, X, and m are 

considered as constants for a given temperature and can be determined 
experimentally fi*om conventional constant stress or load uniaxial tensile 
or compressive creep tests perfonned at constant tenperatures . The cor- 
responding relationship between stress, strain, and time for such tests 
can be found from equation (l) i^ the strain at t = 0 is taken 

as ffo/®! (see fig. 2(a)). Thus 

+ 1 ( 2 ) 

This equation is shown plotted in figure 2(b). 


Adaptation of Stress -Strain -Time Law to 
Bending and Buckling Problems 

When a beam whose cross section can be represented by the idealized 
H-section, shown in figure 3> is bent under the action of combined lateral 
and axial loads, the stresses in each of the flanges may be compressive, 
or the stress in one flange may be conpressive and in the other tensile. 
Thus, in applying equation (l) to the bending problem, the following forms 
of this equation have been used: 


e = (cr/fe;Lc) (3) 

for a flange \inder the action of compressive stress, and hence e and a 
are positive in compression. The subscript c denotes parameters deter- 
mined from conventioneQ. compressive creep tests. Similarly, 

i - -(“Alt) - (■^"‘At) 

for a flange under the action of tensile stresses, and hence e is posi- 
tive in compression and a is positive in tension. The subscript t 
refers to parameters determined from tensile creep tests. 
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Equation ( 3 ) is valid for 0^0 and positive in compression, 
vhereas equation (4) applies vhen 0^0 and positive in tension. 


Stress -Load and Strain-Load Eelations 


If 0c and 0 ^, respectively, are the stresses on the concave and 

convex side of an idealized H-heam loaded in the plane containing the 
web (see fig. 3) and if a moment M and an axially compressive load P 
are applied as indicated in figures 3 and 4, then 


0 ^ = (P/A) + (2M/Ah) 






0 ^ = (p/a) - (2M/Ah) 


(5) 


in which 0 ^ and 0 ^ are considered positive in compression, and 


0 ^ = (2M/Ah) + (P/A) 
0 ^ = (2M/Ah) - (p/a) 


( 6 ) 


in which 0c is positive in compression, while 0 ^ is positive in ten- 
sion. In equations (5) and (6) a/ 2 is the area of each flange and h 
is the distance between flanges . _ 


When both flanges are in compression, the strain in each flange can 
be related to the applied loads with the aid of equations ( 3 ) and (5)« 
Thus ^ 


6^ = (l/E3_^)(3/at)[(P/A) + (2M/Ah)] 
^ = ( l/Eic )( S/St )[( p/a) - (2M/Ah)j 


+ (l/Xc)[(P/A) + (2M/Ah)]“'' 
+ (l/X^)[(P/A) - (2M/Ah)]“‘^ 


(7) 


in which and €^, respectively, are the strain rates on the concave 

and the convex sides of the beam, positive in compression. 
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Similarly, when one flange is in compression and the other in 
tension 

= (l/Elc)(VSt)[(2M/Ah) + (P/A)] + (l/\j[(2M/Ah) + (p/A^ 

k = -(V%t)(^M)[C2M/Ah) - Cp/A)] - (l/>^)[C2M/Ah) - Cp/A)] 

Thus equations (7) are applicable when P/A ^ 2M/Ah and equa- 
tions (8), when P/A ^ 2M/Ah. 


(8) 


Differential Equations of Bending of a Beam Column 

Since the strains and considered positive in compression, 

are related to the radius of curvature p of the centroidal axis of the 
beam by the equation 

(^c - et)/h = l/p ( 9 ) 

the differential equations of bending are readily obtained. Thus, if 
P/A^ 2M/Ah, that is, if both flanges are in compression, equations (7) 
and ( 9 ), together with the condition that for small deflections p is 
related to the deflection w by the relation 

1 / p = ( 9 a ) 

where x is the axial coordinate of a point on the beam axis nnd a sub- 
script X Indicates differentiation with respect to x, yield 

-h( 8 /dt)(wxx) = ( 2 M/Ah^™'^ - 

[(P/A) - (2M/Ahg“^| ( 10 ) 

If p/a ^ 2M/Ah, corresponding to one flange in compression and the 
other in tension, then from equations ( 8 ) and ( 9 ) 

-h(8/8t)(w3^) = (l/Eic)(V^)[(2M/Ah) + (P/A)l + (l/Ei^)(8/8t)[(2M/Ah) - 

(p/a)] + (l/lj[(2M/Ah) + (p/a]“^ + 

(1/XJ[(2M/Ah) - (P/A)]“^ 


( 11 ) 
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Furthermore, if the parameters have the same values in tension and 
compression, the subscript c can he dropped in equation (lO), and 
equation (ll) becomes 

-h(8/8t)(v^) = (yE^Ah)(8M/^) + (l/X) |[(2M/Ah) + (p/A)| + 

[(2M/Ah) - (P/A)[“| (11a) 

Thus, equations (lO) and ( ll) are the differential equations appli- 
cable to the analysis of the creep-deflection behavior of an Idealized 
H-sectlon beam column, the material of which behaves according to equa- 
tions ( 3 ) and (4). 


APPLICATIONS OF DIFFERENTIAL EQUATIONS 
OF BENDING 

Beam Under Pure Bending 


Equation (ll) can be readily applied to- the problem of the determi- 
nation of the creep deflections of an idealized H-section beam under the 
action of constant end couples Mq. Under such conditions this equation 
reduces to 

( 12 ) 

The related boundary and Initial conditions for a beam of length L 
with the origin of the axial coordinate x at one end of the beam are 

that, at X = 0 and L, -w = 0 and, at t = 0, w is 

''o = [(V%c) + CV®it)](«oi'V''i)(*A)[i - (==A)] (13) 

which is obtained from a simple "elastic" analysis. The moment of 
inertia I of the cross section of the beam with respect to the center 

line normal to the web (see fig. 3) is equal to Ah^/4. 

Hence, the solution of equation ( 12) can be expressed as 

^ - ^0 = [(V^c)(2Mo/Ah)“'= + 

(l/^)(2MQ/Ah)”^ (L2/2h)(x/L)[l - (x/L^t (l4) 
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If the material properties of the beam are the same for compressioxi and 
tension, equation (l4) can be expressed as 


w/wq = (E3y\)(2Mo/Ah)™“^t + 1 Cl^a) 


or 

in vhich Oq is the absolute value of the constant stress 2t^/Ah acting 
in each flange. It may be noted that equation (l^b) can also be obtained 
from solution of equations (2), C9)^ aad (9&)* Because the flange 
stresses are constant and equal in magnitude, the ratio of the deflection 
at any time t and the deflection at t = 0 is the same as the ratio 
of the strain at any time t and the strain at t = 0 of a bar under 
the action of constant tensile or compressive stresses (see fig. 2(b) 
and compare eqs. (2) and (l4b)). 


Column With Initial Curvature 

The deflection-time characteristics of an H-section column, whose 
centroidal axis initially deviates from a strai^t line and whose 
defining material parameters are the same for tension and compression, 
can be investigated with the aid of equations (lO) and (lla). For a 
simply supported column with a constant axially compressive end load P, 
the moment M is simply 


M = P(w + Wj^^ 


(15) 


in which w is the deflection due to loads and Wj^ represents the 

initial deviation from straightness Csee fig. 5). Thus, for P/A ^2M/Ah 
and hence, from equation (l5)j for w + ^ h/2, equation (lO) becomes 
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in which 


p/a = a 

w/h = F 

"lA - ^1 ■ 

Equation (l6) is applicable to the present problem so long as both 
flanges are in compression. 

Similarly, for w + wj ^ b/2, that is, for the flange on the concave 
side in compression and that on the convex side in tension, equation (lla) 
yields 

^xx + (oA/e^1)¥ + [(2a)“Ah^] | [f + Fi + (1/2)']“' + 

[f + Fi - (1/2)] “I = 0 (17) 

It may be noted that equations (l6) and (17) are identical if m is 
an odd integer. 

Since the deflection due to loads consists of the time -independent 
deflection Wq, obtained at t = 0, and the time-dependent deflection 

v^, accrued for t > 0, F can be expressed as 

F » F(, + Fq 

Hence , 

F + F^ = F(. + Fij,^ (19) 

in which Fm represents the total deviation of the axis of the loaded 
■‘■o 

column from the x-axis at t = 0 and F^ = w^/h. Hence, equations (l6) 
and ( 17 ) become for F^ + ^ l/2 

Fc^ + (oA/¥jJ)¥c + [(2a) |[(l/2) + F^ + F^ J “ - 

[(1/2) - Fe - Ft^]“| = 0 


( 20 ) 
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and for F_ + Fm ^ I /2 


+ 

^xx 


(aA/Eil)Fc + [(20)“Ah^ j^ljc + Fr^ + ( l/2 j] “ + 


[^c + 


m 


r = 0 


( 21 ) 


If m is iinity, then (since eq. (l) reduces to the stress -strain- 
time law for a Maxwell element) equations (20) and (21) are applicable 
to the analysis of the behavior of the corresponding linearly visco- 
elastic column (see refs. 6, 8, and 9)- If tlie ratio of the initial 
deviations from straightness and the distance between flanges is taken 
as 


Fj^ = fi sin (rtx/l) (22) 

in which L is the length of the beam and the coordinate x is measured 
f 2 x»m a support (see fig. 5)j the deflections of the column can be 
expressed as 


- 1 (23) 

where is the anplitude of the deflection F(, accrued for t > 0, 

fm is the amplitude of the total deviation from the x-axis of the 
■^o 

loaded column at t = 0, and 





The amplitude of the total nondimens ional midspan deflection meas- 
ured from the x-axls is 


% - 


f c + fm 


Equation ( 23 ) shows that for a < Cg the creep deflections of the line- 
arly viscoelastic column become infinitely large only for correspondingly 
large values of time (refs. 8 and 9 )- 
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For values of m other than unity, the collocation method can be 
applied to the solution of equations (20) and (2l), This method will 
be used to investigate the midspan deflection_of the column under the 
assumption that the shape of the original d evi ation from straightness F^ 
is maintained during the bending process. Hence, if F^ is as defined 
in equation (22), then Fg is assumed to be given by the equation 


Fg = f^ sin (nx/L)_ (24) 

Corresponding to equation (22), the total deviation of the axis of the 
loaded column from the x-axis at t = 0 is obtained from an "elastic" 
analysis as 


F™ = f„ sin (kx/l) ( 25) 

■^o -^o 


in which 

Substitution of F,, aaid ^Tq^ respectively, from equations (24) and 

(25) into equations (20) and (2l) and evaluation of the resulting expres- 
sions at X = l/ 2 yield the following equations for the nondimens ional 
midspan deflection f^: For f^, + f^ < l/2 .. 


df, 


,/dr - |[(l/2) + fe + - [(1/2) fc - = 


= 0 (26) 


and for f^ + fip = l/2 


df. 


/dT - fj, + fij^ + (1/2^ + Jfj. + f^^ - (1/2^ ^ = 0 (27) 


in which the nondimens ional time variable t is defined as 


= |Ei(2a)yj4x(aE - a)j|t 


( 28 ) 
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Equation ( 26 ) is applicable if both flanges of the column are in compres- 
sion, whereas equation ( 27 ) applies when one of the flanges is in ten- 
sion. The solution of equation (26) can be written as 


T 


<azy|^[(i/2) + - [(1/2) 



(29) 


in which z = fc + f^^, and hence fji^ ^ z ^ l/2. Thus equation (29) is 
valid provided 0 ^ t < t where corresponds to z = l /2 and hence 

to that time at which one flange is completely unstressed. If f^ ^ l/2, 
then equation ( 29 ) must be applied until t = whereupon the solution 
to equation (27) becomes applicable. This solution is 

^1 + (1/23“"+ - (1/2)]^ (30) 

for 1/2 < z < «> and ^ t ^ where is defined as the criti- 

cal time and corresponds to infinite deflections. Thus, for f^.^ ^ l/2, 

"^cr = -Tl + + (1/2)]“ + t - (1/2)]“| (31) 

and from equation (29) 

^ '"1 = ^ dZy/|[(l/ 2 ) + z]”" - gl/ 2 ) - z]“| ( 32 ) 

If fT^ > 1 / 2 , only one flange is "in con^iression for all values of 

and hence equation ( 27 ) applies throughout the duration of loading. 
Therefore, for f^ > 1/2 


T = 


fj + (1/2)]"+ 6 - (1/2)J"| 


(33) 
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and 0 ^ T ^ Tcr> where 


Ter 



[z - (1/2)] 



( 34 ) 


The integrals in equations (29) to (34) can be evaluated readily for 
integral values of If m is an even integer (or fractional) all six 

equations are required for a complete range of f^i^ and t. If m is 

an odd integer, since equations (26) and (27) are then identical, the 

evaluation of the six equations reduces to the evaluation of equations (33) 

and (34) for all values of fm and t. Thus t-i no longer is of any 

-‘-o 

significance. A summary of equations (29) to (34) is given in table 1 
and the corresponding results of integration for m = 1, 2, 4, and 5 

are given in table 2. The integrations indicated in equations (29) to 
(34) are performed in reference 11 in terms’ of m, vhere m is any 
integer. It may be noted from comparison of equation (23) with the 
results in table 2 that, for m = 1, the collocation method yields the 
exact results of equation (23). 


The expressions in table 2 were used to obtain the deflection-time 

curves shown in figures 6(a) to 6(d). These curves relate the nondimen- 

sional time -dependent midspan deflection f^ (accrued for t > o) to 

the time variable t and the initial total deviation from the x-axis of 
the midspan of the loaded column figure 6(a) the curve designated 

m = 2, 0 g T ^ T]_> or m = 1, 0 ^ t ^ w, is valid for the case in which 

m = 2 and both flanges are in compression, as well as for the entire time 

range for m = 1 (see table 2). In table 3 the time to failure, charac- 
terized by T^j., and the time to zero stress in the "tension" flange, 

represented by T;,, are listed for the m and fm values considered 

in figures 6(a) to 6(d). 

The curves in figures 6(a) to 6fd) show that for the same value of 
the exponent m > 1 the creep deflections becon» large much more rapidly 
for large values of frp^ than for small values of this parameter. The 

result is, of course, that the critical time parameter decreases 

as f^i^ increases (see table 3)- These results are as expected, since 
large values of frp^ indicate either large initial deviations of the 

column from straightness or large values of the average axial compressive 
stress a < Og. From table 3^ as well as from a comparison of the loca- 
tions of the asymptotes shown in figures 6(a) to 6(d), it is seen that 
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for large values of the critical time parameter decreases 

vith increasing m, whereas for small values of fm this parameter 

■*■0 

increases with m. In this respect it may he noted, however, that since 
is proportional to it is not a true measure of the actual 

critical time t^r (eq. (28)). Thus, it is probable that for real, 
materials the actixal critical time t^.^ decreases monotonlcally with 
increasing m, since t is inversely proportional to b^. 

DISCUSSION 


Differential equations of bending of an Idealized H-sectlon beam 
column were derived for a non linearly viscoelastic material whose mechan- 
ical behavior is analogous to a model consisting of a spring connected in 
series with a dashpot. The spring was assumed to be elastic with a modu- 
lus corresponding to the strain intercept of the projection of the second- 
ary region of a conventional tensile or compressive creep curve, whereas 
the dashpot was considered as nonlinear with a strain rate proportional 
to a power of the stress. Thus, the conventional creep curve is repre- 
sented by a strai^t line, the slope of which is the secondary creep rate 
of the actual creep curve. The general equations, one which applies when 
both flanges are in compression and the other when one flange is in com- 
pression and the other in tension, are applicable to materials whose 
mechanical properties for tension differ from those for compression. 

The results of the einalysis of the creep deflections of a beam under 
pure bending showed that the deflections vary linearly with time. This 
is to be expected, since the stress in each flange remains constant with 
time, and hence each flange behaves in a manner analogous to a bar under 
constant stress . If the material properties of the beam are the same for 
tension and compression, the ratio of the deflection at any time t to 
the deflection at t = 0 is identical to the ratio of the strain at the 
same time t to the strain at t = 0 of a bar under the action of a con- 
stant stress whose magnitude is that of tlte stress in either flange of the 
beam (see eqs. ( 2 ) and (l4b) and fig. 2(b)). 

The nonlinear differential equations for the bending of a simply 
supported H-section column with Initial sinusoidal deviation from strai^t- 
ness were solved with the aid of the assumption that the shape of the 
loaded column remained sinusoidal. The material properties of the beam 
were assumed to be the same in tension and compression. In contrast with 
the results presented in references 8 and 9 for linearly viscoelastic 
columns, the nonlinear ly viscoelastic columns investigated in the present 
paper possessed "critical times." Thus for stresses less than the Euler 
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stress the deflections of an initially crooked column become infinitely 
large in a finite time. In figures 6 (a) to 6 (d) curves are presented 
showing the effect of the initial deviation from straightness of the loaded 
column on the time behavior of the additional deflections obtained with 
time. As expected, for a given exponent in tte power law for viscosity the 
creep deflections increased more rapidly with time for the larger values of 
amplitude of the initial deviation from straightness than for the smaller 
values. Thus, the critical time decreases very rapidly with increasing 
initial deviation from straightness and end load. 


of Brooklyn, 

July 23, 1952. 


Polytechnic Institute 
Brooklyn, N. Y., 
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TABI£ 2.- SOLCWIOH OF CHBEP-BUCKLIHG EQUATIOHS FOR m « 1, 2, 3j AHD 5 - Ccncluded 
(b) m odd integer ^0 < < a>, ^ z oaj 


Solution 

’’’cr j 

a «= 1, t/r = 

^ » % e^" 
■*•0 

00 

m = 3, t/r - X(oe - ^)/(2E]^o3) 

f r 

- e3T) + 3 

i log / 1 + — 2 — \ 

' 1 v) 

m = 5, t/r » x(og - 

T ^ log -5 /£ + log + ? - ^^^5 _ / 2 _ 1 \ log + ? + 2|/^_ 

5 %o '5 ‘'5/ l^fTo^ + 5 - 2)/5 ^5 ^ / 4f^^2 + 5+2^ 

(1 * w) ('‘*' 0 " * 5 - ^) t 

(I ■ ^) (V - 5 • =1^) - 1 >« 


• t 
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TABLE 3.- VALUES OF AND CORRESPOKDING TO VALUES 

OF m. AND OF FIGUEES 6Ca) TO 6Cd) 


m = 2 

m = 3 

% 

■^o 

•"l 

T 

cr 

■"cr 

0.01 

1.956 

2.74i 

2.974 

.05 

1.152 

1.937 

1.902 

.10 

.8045 

1.590 

1.444 

.20 

.4582 

1.244 

.9943 

.50 

0 

.7854 

.4621 

.80 


.5586 

.2586 

1.00 


.4637 

.1865 

2.00 


.2450 

.05728 

m = 4 

m = 5 

^T 

-■-o 

^1 

^cr 

^cr 

0.01 

3.565 

3.891 

5.612 

.10 

1.282 

1.608 

1.985 

1.00 


.09112 

.04767 
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(a) Experimental. 
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(b) Idealized. 

Figure 2.- Experimental and idealized creep curves. 





Figure 6.- Variation of deflections accrued for t > 0 wltli time 
parameter T for m = 2, 5^ sind 5 S'S‘3. several values of 

initial-deflection parameter fy . 



















